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Abstract 

The Tokai-to-Kamioka-and-Korea (T2KK) neutrino oscillation experiment under 
examination can have a high sensitivity to determine the neutrino mass hierarchy for 
a combination of relatively large (~ 3.0°) off-axis angle beam at Super-Kamiokande 
(SK) and small (~ 0.5°) off-axis angle at L ~ 1,000 km in Korea. We elaborate 
previous studies by taking into account smearing of reconstructed neutrino energy due 
to finite resolution of electron or muon energies, nuclear Fermi motion and resonance 
production, as well as the neutral current vr'' production background to the i'^ — > 
oscillation signal. It is found that the mass hierarchy pattern can still be determined 
at 3cj level if sin^ 26Irct = 4|C/e3p(l - |t/e3p) ^0.08 (0.09) when the hierarchy is 
normal (inverted) with 5 x 10^^ POT exposure, or 5 years of the T2K experiment, if a 
100 kton water Cerenkov detector is placed in Korea. The vr'^ backgrounds deteriorate 
the capability of the mass hierarchy determination, whereas the events from nuclear 
resonance productions contribute positively to the hierarchy discrimination power. 
We also find that the vr'^ backgrounds seriously affect the CP phase measurement. 
Although 5mns can still be constrained with an accuracy of ~ ±45° (±60°) at la level 
for the normal (inverted) hierarchy with the above exposure if sin^ 20rct ^ 0.04, 
CP violation can no longer be established at 2>a level even for Jmns = ±90° and 
sin^ 26'rct = 0.1. About four times higher exposure will be needed to measure 5mns 
with ±30° accuracy. 



*E-mail address : naotoshi(S)post. kek.jp 



1 Introduction 



The SNO experiment found that the from the sun changes into the other active neutri- 
nos [1]. The atmospheric neutrino observation at SK reported that z/^ and z/^ oscillate into 
the other active neutrinos [2]. Recently, the MiniBooNE experiment [3] reported that the 
LSND [4] observation of rapid z/^ oscillation has not been confirmed. Consequently, 

the three active neutrinos are sufficient to describe all the observed neutrino oscillation 
phenomena. 

Under the three generation framework, neutrino fiavor oscillation [5, 6] is governed by 
2 mass-squared differences and 4 independent parameters in the MNS (Maki-Nakagawa- 
Sakata) matrix [5], that is 3 mixing angles and 1 CP phase ((5mns)- The absolute value 
of the larger mass-squared difference, l^m^g] and one combination of the MNS matrix 
elements sin^2^ATM = 4|t/^3p(l — |t/^3p), are determined by the atmospheric neutrino 
observation [7, 8, 2, 9], which have been confirmed by the accelerator based long base- 
line neutrino oscillation experiments K2K [10] and MINOS [11]. However, the sign of 
^mfg = m| — has not been determined. Both the magnitude and the sign of the smaller 
mass-squared difference 5m^2 — ^2 ~ "^i; another combination of the MNS matrix ele- 
ments sin^2^soL = 4|C/eiC/e2p are determined by the solar neutrino observations [12, 1] 
and the KamLAND experiment [13]. The last independent mixing angle (^rct) has not 
been measured yet, and the reactor experiments [14] give upper bound on the combination 
sin^2^j^cT = 4|C/e3p(l — I^e3p)- The leptonic CP phase, (5mns — — argC/gs [15], is unknown. 

There are many experiments which plan to measure the unknown parameters of the three 
neutrino model. In the coming reactor experiments. Double CHOOZ [16], Daya Bay [17], 
and RENO [18] plan to measure the unknown element |C/e3| from the Pg survival probability. 
The Tokai-to-Kamioka (T2K) neutrino oscillation experiment [19], which is one of the next 
generation accelerator based long baseline experiments, also plans to measure jC/gsl by 
observing the v^^ — > v^, transition event, whose rate is proportional to ^ezU^z^ ■ 

However, the sign of ^vn\-^. or the mass hierarchy pattern, will remain undetermined even 
after these experiments. It is not only one of the most important parameters in particle 
physics but also has serious implications in astronomy and cosmology. For instance, if ^m^g 
is negative (inverted hierarchy), the prospects of observing the neutrino-less double beta 
decay are good, while the matrix element \Ue'i\ is affected by quantum corrections such 
that its high energy scale value depends on the Majorana phases [20] in the large tan/? 
supersymmetric See-Saw scenario [21]. In astronomy, the mass hierarchy pattern affects 
the light elements synthesis in the supernova through neutrino-nucleon interactions; the 
yields of ^Li and ^^Be increase for the normal hierarchy (^^13 > 0) if sin^ 2^rct ^ 10~^ [22]. 
In cosmology, the dark matter content of the universe depends on the mass hierarchy. 

In the previous studies [23, 24, 25], we explored the physics impacts of the idea [26] of 
placing an additional far detector in Korea along the T2K neutrino beam line, which is now 
called as the T2KK (Tokai-to-Kamioka-and-Korea) experiment. In particular, we studied 
semi-quantitatively the physics impacts of placing a 100 kton water Cerenkov detector 



in Korea, about 1000 km away from J-PARC (Japan Proton Accelerator Research Com- 
plex) [27], during the T2K experiment period [19], which plans to accumulate 5 x 10^^ POT 
(protons on target) in 5 years. We find that the neutrino-mass hierarchy pattern can be 
determined by comparing the z/^ z/g transition probability measured at SK (L = 295 km) 
and that at a far detector in Korea [23], if sin^ ^Oj^^t ^ 0.05 for 3a. The CP phase can also 
be measured if sin^ 29^ct ^ 0.02 with ±30° accuracy, since the amplitude and the oscillation 
phase of the transition probability arc sensitive to sin (^MNg and cos^nj^g, respec- 

tively [23, 24]. We also find that the octant degeneracy between sin^ ^atm = 0.4 and 0.6 for 
sin^2^ATM = 0.96 can be resolved if sin^ 2^rct ^0.12 [25]. In the above studies [23, 24, 25], 
a combination of 3.0° OAB (off-axis beam) at SK and 0.5° OAB at L = 1000 km in Korea is 
found to be most efficient, mainly because of the hard neutrino spectrum of the 0.5° OAB. 
In alternative studies [28] of the T2KK setup, an idea of placing two identical detectors 
at the same off-axis angle in Kamioka and Korea has been examined. The idea of placing 
far and very far detectors along one neutrino baseline has also been studied for the Fermi 
Lab. neutrino beam [29]. 

The T2KK experiment has a potential of becoming the most economical experiment 
to determine the mass hierarchy and the CP phase, if sin^ 2^rct is not too small. In 
this paper, we re-evaluate the T2KK physics potential by taking into account smearing 
of the reconstructed neutrino energy due to finite resolution of electron or muon energies 
and the Fermi motion of the target nucleon, as well as those events from the nuclear 
resonance production which cannot be distinguished from the quasi-elastic events by water 
Cerenkov detectors. We also study contribution from the neutral current 7r° production 
processes which can mimic the appearance signal. 

This article is organized as follows. In section 2, we fix our notation and give approx- 
imate analytic expressions for the neutrino oscillation probabilities including the matter 
effect. The relations between the experimental observables and the three neutrino model 
parameters are then explained by using the analytic formulas. In section 3, we show how 
we estimate the event numbers from the charged current (CC) and neutral current (NC) 
interactions by using the event generator nuance [30]. In section 4, we present the 
function which we adopt in estimating the statistical sensitivity of the T2KK experiment 
on the neutrino oscillation parameters. In section 5, we show our results on the mass 
hierarchy determination. In section 6, we show our results on the CP phase measurement. 
In section 7, we give the summary and conclusion. In Appendix Appendix A, we present 
a parameterization of the reconstructed neutrino energy distribution as a function of the 
initial neutrino energy for CCQE and resonance events. 

2 Notation and approximate formulas 

In this section, we fix our notation and present an analytic approximation for the neutrino 
oscillation probabilities that is useful for understanding the physics potential of the T2KK 



experiment qualitatively. 



2.1 Notation 

The neutrino flavor eigenstate |z/a) {a = e,fj.,T) is a mixture of the mass eigenstates |z/j) 
{i — 1,2,3) with the mass as 

3 

Wa) ^^Uai Wi) , (1) 

i=l 

where U is the unitary MNS (Maki-Nakagawa-Sakata) [5] matrix. We adopt a convention 
where Uei, C/e2, C/^3, > and 5mns = -argC/e3 [15, 31]. The 4 parameters, Ue2, C/^3, 
lUesl, and 5mns, can then be chosen as the independent parameters of the 3x3 MNS matrix. 
All the other elements are determined uniquely by the unitarity conditions [31]. 

The atmospheric neutrino observation [7, 8, 2, 9] and the accelerator based long base- 
line experiments [10, 11], which measure the survival probability, are sensitive to the 
magnitude of the larger mass-squared difference and sin^2^ATM [H]: 

sin^ 2^ATM > 0.90 (90% C.L.) , (2a) 
\6ml^\ = (2.43 ± 0.13) X lO^W. (2b) 

The reactor experiments, which observe the survival probability of z/g at L ~ 0(1) km, 
are sensitive to If^mfgl and sin^2^jj^cT- The CHOOZ experiment [14] finds 

sin2 20RCT < (0.20, 0.16, 0.14) (3a) 
for ISmj^l = (2.0, 2.5, 3.0) X 10-^eV^ (3b) 

at the 90% confidence level. 

The solar neutrino observations [12], and the KamLAND experiment [13], which measure 
the survival probability of and u^, respectively, at much longer distances are sensitive to 
the smaller mass-squared difference, 5m^2) ^e2- The combined results [13] find 

sin2 2^soL = 0.87 ±0.04, (4a) 
5ml^ = (7.59 ± 0.21) x lO'^eV^ . (4b) 

The sign of 5m^2 has been determined by the matter effect inside the sun [32]. 

With a good approximation [33], we can relate the above three mixing factors, eqs. (2a), 
(3a), (4a) with the elements of the 3x3 MNS matrix; 

sin 6'atm = Ui^3 = sin $23 cos 9i3 , (5a) 
sin6lRCT = |C^e3l = sin6li3 , (5b) 
sin26'soL = 2C/eiC/e2 = sin26li2COS^6'i3, (5c) 

where the three mixing angles Oij = 9ji are defined in the region < 612, ^13, ^23 < 7r/2 [15]. 
In the following, we adopt 

sin^ATM, sin^jjCT, and sin^^soL defined above as the indepen- 
dent real mixing parameters of the 3x3 MNS matrix. 



2.2 Approximate formulas 



The probability that an initial flavor eigenstate with energy E^, is observed as a flavor 
eigenstate |z/^) after traveling a distance L in the matter of density p{x) (0 < a; < L) along 
the basehne is 



(i//3|exp ^— i j H{x)dx^ 



(6) 



where the Hamiltonian inside the matter is 
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Here Gp is the Fermi constant, E^y is the neutrino energy, ne{x) is the electron number 
density, and p{x) is the matter density along the baseline. In the translation from ne{x) 
to p{x), we assume that the number of the neutron is same as that of proton. To a good 
approximation [19, 34], the matter proflle along the T2K and T2KK basehnes can be 
replaced by a constant, p{x) = po, and the probability eq. (6) can be expressed compactly 
by using the eigenvalues (A^) and the unitary matrix If of eq. (7); 

Pu.-^u, = 6^fs- A J2MK^Up^u,Ju;^) sin' ^ + 2 J2Mu:^^^^ 
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All our numerical results are based on the above solution eq. (9a), leaving discussions of 
the matter density profile along the baselines to a separate report [34]. Our main results 
are not affected significantly by the matter density profile [34] as long as the mean matter 
density along the baseline (po) is chosen appropriately. 

Although the expression eq. (9a) is not particularly illuminating, we find the following 
approximations [23, 24] useful for the T2KK experiment. Since the matter effect is small 
at sub GeV to a few GeV region for p ~ 3 g/cm^, and the phase factor A12 in the vacuum, 
where 

2 2 
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is also small near the first oscillation maximum, |Ai3| ~ tt, the approximation of keeping 



the first and second order corrections in the matter effect and A12 [35, 23, 36, 24] 



(11a) 



4 sin^ ^ATM sin^ O^ct <! (1 + sin^ ( ^ ) + 5^ sin A13 } + C% (lib) 



has been examined in ref. [24]. Here A^^ and B^^ are the corrections to the amphtude and 
the osciUation phase, respectively, of the survival probability. When 1^"^! and are 
small, eq. (lib) reduces to 

P,^^,, ^ 4 sin^ sin^ ^^ct (1 + A') sin^ + B'^ + C , (12) 

similar to the survival probability, eq. (11a). We therefore refer to B^ in eq. (lib) as 
the oscillation phase-shift, even thought it can be rather large (~ 0.4). 

For the survival probability, eq. (11a), it is sufficient to keep only the linear terms in 
A12 and a, 

A- - --^^^^sin^^^ex, (13a) 
A^^E cos2 ^ATM 

-B^ = ^ ^"'^^ ATM ^.^2 0^^^ 12 ^^^g2 0^^^ _^ ^^^2 0^^^^ gjj^2 Q^^^ Q^^^ 

4:11/ COS l/ATM ^ 

-tan6'ATMsin26'soLsin6'RCTCOs5MNs) ■ (13b) 

The above simple analytic expressions reproduce the survival probability with 1% accuracy 
throughout the parameter range explored in this analysis, except where the probability is 
very small, {P^^^^^ ^10~^). In eq. (13a), the magnitude of A^^ is much smaller than 
the unity because of the constraints (2a) and (3a), and hence the amplitude of the 
survival probability is not affected significantly by the matter effect. This means that 
sin^ 2^ATM can be fixed by the i/^ disappearance probability independent of the neutrino 
mass hierarchy and the other unconstrained parameters. The phase-shift term B^ affects 
the measurement of l^m^gj. However, the magnitude of this term is also much smaller 
than that of the leading term, A13/2, around the oscillation maximum IA13I ~ tt, because 
cos26'atm = -x/l - sin^ 26'atm < by eq. (2a) and A12/A13 < 1/30 by eqs. (2b) and 

(4b). The smallness of the phase shift term does not allow us to determine the sign 
of Ai3, or the neutrino mass hierarchy pattern, from the measurements of the survival 
probability only. 

For the i/^ — > Ug transition, eq. (lib), we need to retain both linear and quadratic terms 



of Ai2 and a to obtain a good approximation; 
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+ f ^ J sin^ 6'rct sin^ 6'atm ■ ( 14c) 

Here, the first and second terms in eqs. (14a) and (14b) are the hnear terms of A12 and a 
respectively, while the other terms and all the terms in eq. (14c) are quadratic in A12 and 
a. These quadratic terms can dominate the oscillation probability when sin^ Q^q^. is very 
small. We find that these analytic expressions, eqs. (lib) and (14), are useful throughout 
the parameter range of this analysis, down to sin^ ^rct = 0, except near the oscillation min- 
imum. The amplitude of the i/^ — > Vf. transition probability, is sensitive to the mass 
hierarchy pattern, because the first term of changes sign in eq. (14a), with cos 2^rct ~ 1- 
When L/E is fixed at |Ai3| ~ tt, the difference between the two hierarchy cases grows with 
L, because the matter effect grows with E; see eq. (8). The hierarchy pattern can hence 
be determined by comparing Py^^y^ near the oscillation maximum |Ai3| ~ tt at two vastly 
different baseline lengths [23, 24]. 

Once the sign of A13 is fixed by the term linear in a, the terms linear in A12 allow us 
to constrain sin^MNs via the amplitude A^, and cos^mns via the phase shift B^. Therefore, 
5mns can be measured uniquely once the mass hierarchy pattern and the value of sin^ 2^^jiCT, 
which may be measured at the next generation reactor experiments [16, 17, 18], are known. 



3 Signals and Backgrounds 

In this section, we show how we estimate the event numbers from the charged current (CC) 
and the neutral current (NC) interactions. First, we explain how the signal CCQE events 
are reconstructed by water Cerenkov detectors, and study contributions from the inelastic 
processes when none of the produced particles emit Cerenkov lights and hence cannot be 
distinguished from the CCQE events. Next in subsection 3.2, we study NC production of 
single 7r°, which can mimic the z/^ v^. appearance signal when the two photons from 7r° 
decay cannot be resolved by the detector. Finally, we show the sum of the signal and the 
background events. 




Figure 1: Reconstructed energy distributions for z/^ CC events on the water target at 
Ejj = 1 GeV (a) and = 2 GeV (b) according to the event generator nuance [30] when the 
/i momenta are measured exactly. Among 10^ generated events about 73% are CC events 
at both energies, which consist of CCQE events, nuclear resonance production events, and 
the others including deep-inelastic processes. After the CCQE selection cuts of eq. (16) 
are applied, the blue shaded region survive. 

3.1 CC events 

In accelerator based long baseline experiments, one can reconstruct the incoming neutrino 
energy E^^ by observing the CCQE events {vg^n ip or u^p in) if the charged lepton 
{i = fi or e) momenta are measured and the target nucleons are at rest, since the neutrino 
beam direction is known. In practice, however, the lepton momentum measurements have 
errors, the nucleons in nuclei have Fermi motion, and some non-CCQE events cannot be 
distinguished from the CCQE events. None of those uncertainties has been taken into 
account in the previous studies of refs. [23, 24, 25]. In this and the next subsections, we 
study them for CC and NC processes, respectively, for a water Cerenkov detector by using 
the event generator nuance [30]. 



3.1.1 Event selection 

In a CCQE event, Ufji ip, the neutrino energy E^ can be reconstructed as 

_ rUnEe - mj/2 - {ml - m^) /2 

-^rec — ^ — ; 7 , 

nin — El + pi COS U 

in terms of the lepton energy {Ei), total momentum {pe), and its polar angle 6 about the 
neutrino beam direction, if a target neutron is at rest. For an anti- neutrino CCQE event, 
uip — > in, rup and m„ should be exchanged in eq. (15). 

In reality, the target nucleons inside nuclei has Fermi motion of about 100 MeV, and the 
measured e and /i momenta have errors. Therefore, i^rec of eq. (15) is distributed around 



the true Ei,, even for the CCQE processes. 

The CCQE events are selected as 1-ring events in a water Cerenkov detector by the 
following criteria [10, 19] : 

Only one charged lepton (i = ji^ or e^) with \pi\ > 200 MeV, (16a) 

No high energy vr^ (|p^±| > 200 MeV) , (16b) 

No high energy 7 {\p^\ > 30 MeV) , (16c) 

No 7r°, Ks, Kl, and . (16d) 

The lower limit of the total momentum in the first criterion in eq. (16a) is from the 
threshold of the water Cerenkov detector ior i — n [8]. tt^ with \p\ > 200 MeV or 7 with 
\p\ > 30 MeV gives rise to an additional ring. Also, 7r°, Kl, Ks, and K"^ always decay 
inside the detector, making additional rings. 

Figure 1 shows the E^ec distribution of the CC events at E',^ = 1 GeV (a) and 
El, — 2 GeV (b) on the water target, according to the event generator nuance [30]. Among 
the 10^ events at each energy, about 73% are CC events (the rests are NC events) which 
consist of CCQE events, nuclear resonance production, and the others including deep 
inelastic events. After the CCQE selection cuts of eq. (16) are applied, the blue shaded 
region survives, which consists of the CCQE events and the other events where the produced 
77=*= are soft. We call the non-CCQE events which survives the selection cuts of eq. (16) 
"resonance events", since most of them come from single soft tt^ emission from the A 
resonance. The CCQE events and the resonance events are observed as two peaks in 
the reconstructed energy which are separated by about 380 MeV at E^,^ ~ 1 GeV, rather 
independent of the initial energy. This is because the origin of the distance between 
the two peaks mainly comes from the mass difference between the nucleon and the A 
resonance, which scales as 

^J<Z<L ^ 340 MeV (17) 

2m„ ^ ' 

in eq. (15). Because the peak value of the factor, — p^cosO, in the denominator of 
eq. (17) decreases from about 100 MeV at E^ = 1 GeV to about 50 MeV at E^, = 2 GeV, 
the difference in the peak locations decreases slightly from about 380 MeV at E^, = 1 GeV 
in Fig. 1(a) to about 360 MeV at = 2 GeV in Fig. 1(b). The half width of the CCQE 
peak is about 60 MeV, almost independent of E,^, because it comes from the Fermi motion 
of the target nucleons inside nuclei. 



3.1.2 Lepton momentum resolutions 

After selecting the CCQE-like events, we examine the detector resolution which further 
smears the E^-^c distribution. Wc use the momentum and angular resolutions of the muon 
and electron at SK [8], which are shown in Table 1. For the momenta around 1 GeV, 
the momentum resolutions are about a few % and the angular resolutions are about a few 
degrees for both /x and e. 
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Table 1: The momentum and angular resolution of /i- and e-momenta at SK [8]. 
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Figure 2: ii^j-ec distribution of the CC events on the water target for monochromatic energy 
z/^ at Eu = 1 GeV (a) and at E,^ = 2 GeV (b), after the CCQE selecting cuts, generated by 
nuance [30]. The dotted curves show the distributions when the /i^ momenta are measured 
exactly, the boundaries of the blue region in Fig. 1, whereas the solid lines show the 
distributions after the finite momentum resolution of Table 1 is taken into account. 

In Fig. 2, we show by solid curves the Ej-ec distributions after taking account of the fi^ 
momentum resolutions of Table 1, while the dotted lines show the distributions when the 
/i^ momenta are measured exactly, which are the boundaries of the blue shaded region in 
Fig. 1. 

The total width of the CCQE peak is now the sum of the effects from the Fermi motion 
(cpcrmi), the momentum resolution (asp/p), and the angular resolution (ase); it grows with 
E^, because (Xsp/p grows with the lepton momentum. For instance, the half width is about 
60 MeV for Ey = 1 GeV and 70 MeV for Ey = 2 GeV. As a consequence of the energy 
dependence for the total width, the peak height of the CCQE events becomes lower, by 
about 80% for E^ = 1.0 GeV and 67% for 2.0 GeV. 

The ii^rec distribution for the z/g CCQE events are very similar, and we do not show them 
separately. Small differences, due to poorer momentum resolution of electrons in Table 1, 
are reflected in our parameterizations in the next subsection. 



3.1.3 Pctrameterization for the CCQE events 

In this section, we present our parameterization of the E^ec distribution of the CCQE events 
for a given initial or i/g energy E^,, after taking account of the /i- and e- momentum 
resolutions of Table 1. 

The Erec distribution from the CCQE events can be reproduced by three Gaussians, 

/-^-(^^; - i4;y g (- t^ET'^'l ' '''' 

where the index a is for fj, or e, with r^{Ei,) — 1. The factor A°'{E^) ensures the normal- 
ization 

J /„^°Q^(£;rec;^.)^i^rec = l. (19) 

The variance a", the energy shift (n = 1,2,3), and the coefficients, rf" and rg, are 
parameterized as functions of the incoming neutrino energy E,^. These parameters depend 
on the neutrino species, or z/g, because of the mass difTcrcncc in cq. (15), the difference in 
the momentum resolutions in Table 1, and also because of small differences in the CC cross 
sections at low energies [30]. Our parameterization* is given in Appendix A.l, cqs. (A5)- 
(AlO) which is vahd in the region 0.3 GeV < E^ < 6.0 GeV and 0.4 GcV < E,^^ < 5.0 GcV 
for both z/^ and Ug. For the sake of keeping the consistency with the previous studies in 
ref. [23, 24, 25], those events with E^ec < 0.4 GeV are not used in the present analyses. 

In Fig. 3, we show the E^ec distribution of the CCQE events. The solid circles show the 
distributions generated by nuance [30], and the histograms show our smearing functions of 
eq. (18). Figures. 3(a) and (b) are for and fe, respectively, E^ — 1 GeV, and (c) and 
(d) are for those ai E^, — 2 GeV. The area under each distribution is normahzed to unity. 

3.1.4 Nucleeir resonance contributions 

The E^cc distribution of the non-CCQE events which pass the CCQE selection cuts of 
eq. (16) can also be parameterized. Most of them come from the A resonance production, 
and the resonance peak in the E^-^c distribution is observed in Figs. 1 and 2. For E,^ < 
1.2 GeV, 3 Gaussians suffice to reproduce the E'rec distributions generated by nuance [30]; 

/r(^„; < 1.2GeV) = t KiE.) exp (- . (20) 

while at high energies {E^, > 1.2 GeV), we need 4 Gaussians, because the number of 
contributing resonances grows with Ei,; 

/r(^,„; E. > 1.2GeV) ^ ^ P'JE.) exp . (21) 

*A computer code (C/C++) for the parameterization are available from the authors, or directly from 
the web-site [37]. 
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Figure 3: Normalized i^rec distribution of the CCQE events (solid circles) and the resonance 
events (solid diamonds) after the CCQE selection cuts of eq. (16) and the momentum 
resolutions for /i and e in Table 1 are applied. The events are generated by nuance [30]. 
The solid line shows our parameterization for the CCQE events and the dotted line is for 
the resonance events: (a) z/^ and (b) at = 1 GeV; (c) and (d) z/g ai E,, = 2 GeV. 



Around E,^ ~ 1.2 GeV, both parameterizations are valid. Here again a is // or e, f^'^{Ei,) — 
f^'^{Ey) = 1, and the factors A{E^) and A{Et,) assure that the smearing functions are nor- 
mahzed to 1 as in eq. (19). The variances 6"" and cr", the energy shifts 6E", 5E", and 
the relative normalization factors and (n 7^ 1) are all parameterized as functions of 
the incoming energy E^, which are given in Appendix A. 2. The shape of the E^^c distri- 
bution for the "resonance" events are also shown in Fig. 3. The solid diamonds show the 
distribution of non-CCQE events generated by nuance [30] after the CCQE selection cuts 
of eq. (16) and the momentum resolutions of Table 1 are applied. The dotted histograms 
show our smearing functions, eqs. (20) and (21). 

3.2 NC events 

The key observation of ref. [23, 24] for the T2KK proposal is that it is advantageous to 
observe the first oscillation maximum ([Ai3[ ~ tt) at two vastly different baseline lengths, 
L — 295 Km at SK and L ~ 1000 km in Korea. Higher energy neutrino beam, or small off- 
axis angle, is hence desired for the far detector in Korea. However, the use of high energy 
(broad band) beam gives rise to a serious background for the — > oscillation signal. 
The single 7r° production via the neutral current (NC), whose cross section grows with E^, 
cannot always be distinguished from the i/^ — > Ve signal in a water Cerenkov detector. In 
this subsection, wc study the NC 7r° production background in detail and estimate its E^^ 
distribution by using the momentum distribution of misidentified tt'^'s. 

3.2.1 Event selection 

We simulate the NC 7r° production background as follows. By using the neutrino flux^ 
of the T2K beam at various off-axis angles between O.0°(on-axis) and 3.0°, and by using 
the total cross section (t^^^. [a — v^, P^, v^, and Pg) off the water target [30], both CC and 
NC events are generated by nuance [30] for a water Cerenkov detector of 100 kton fiducial 
volume at L = 1000 km, with 5 x 10^^ POT. All the generated events are then confronted 



against the following selection criteria : 

No charged leptons. (22a) 

Only one 7r° . (22b) 

No high energy tt^ (|p^±| < 200 MeV.) (22c) 

No high energy 7 (|p^| < 30 MeV.) (22d) 

No Kl, Ks, K^. (22e) 



The first condition, eq. (22a), selects NC events, and the others eliminate multi-ring events. 
The 7r° momentum distribution after the above cuts is shown in Fig. 4(a) for various off- 
axis beams. We find that the number of single 7r° events grows with decreasing off-axis 

t All the on- and off- axis neutrino flux distributions of the T2K beam used in this report are available 
from the authors, or directly from the web-site [37]. 
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Figure 4: (a): The 7r° momentum distribution of the single 7r° NC events selected by the 
criteria, eq. (22), at various off-axis angles. The event numbers are obtained for a 100 kton 
water target at L = 1000 km with 5 x lO^^POT, according to nuance [30]. (b): Probability 
that a 7r° cannot be distinguished from e^, according to eq. (26). The common horizontal 
axis measures the 7r° momentum. 

angle, especially for the angles below 2.0° which have been envisaged in ref. [23, 24, 25] as 
an optimal choice for the far detector in Korea. 



3.2.2 



misidentification probability 



Figure. 4(a) shows that there are many single- 7r° events from the NC interactions, especially 
for smaller off-axis angles. Some of them become backgrounds of the oscillation 
signal, because the two photons from 7r° are not always resolved by a water Cerenkov de- 
tector. When one of the two photons is much softer than the other, the soft photon dose 
not give a clear ring, resulting in a single-ring (e-like) event. In addition, when the photons 
have a small opening angle the overlapping rings cannot always be resolved. 

We therefore parameterize the probability of misidentifying vr'' as an in terms of the 



energy ratio and the opening angle of the two photons in the laboratory frame. The energy 
fraction of the softer photon in the laboratory frame 

x = — ^ {E,<E^) (23) 

can be expressed as 

x-^(l-/5cos^) , (24) 

in terms of the smaller polar angle (cos^ > 0) of the photon momentum in the 7r° rest 
frame about the polar axis along the 7r° velocity (/9) in the laboratory frame. The opening 
angle between the two photons in the laboratory frame is then 

cos^,, = l-2^^. (25) 

It is clear from eqs. (24) and (25) that when the 7r° momentum is relativistic {(3^1) either 
one of the photons becomes soft {x <^ 1) around cos^ ~ 1, or the two photons become 
coUinear, cos 6'^^ ~ 1. 

By using the energy fraction x and cos^-y-y, the n^-e^ misidentification probability can 
be parameterized as 



Pe/7ro(b7ro|) — / [6 (a;° — x) + ©(x — X°)©(cOS — COS ^^"^)/(x, COS ^^^)] d COS ^ , 

Jq 



(26) 



where Q{x) is the step function. The first step function in the r.h.s. tells that the 7r° 
is misidentified as an when the energy fraction x of the soft photon is smaller than 
When both photons are hard {x^ < x < 0.5), it is still misidentified as an when 
cos^^^ > cos^^^. We introduce a fudge factor 

/(X, cos e,,) = 1.0 - f """" 'l ( , (27) 

in order to take account of detector performance. We show in Fig. 4(b) the ii^-e^ misiden- 
tification probabihty, Pe/7ro(b7ro|), of eq. (26) for x^ — 0.2 and — 17°, which reproduces 
quahtatively the typical performance of water Cerenkov detectors. The leadoff energy, 
IPtt"! — 0.1 GeV, and the height of the plateau, Pe/Tr" — 0.4, are dictated by the first step 
function in eq. (26), which tells that the two photons are not resolved when the softer 
photon has an energy fraction less than 0.2. The second term in eq. (26) determines the 
kink structure around \pT^a\ — 0.9 and 11.0 GeV, as well the asymptotic behavior at high 
7r° momentum. 

The number of the e-like events from the tt'^ background can now be calculated as the 
product of the vr'' event number in Fig. 4(a) and the probability Pe/Tr" ( IPtt" I ) in Fig. 4(b). 
The reconstructed energy E'rec of each tt'' background event is calculated from the 7r° energy 
and the scattering angle by assuming the electron mass. 



3.3 The event numbers 



We calculate the numbers of and Ug CC events from the primary and the secondary 
beam in the i-th energy bin, E^^^ < E < El^^, as 



u^)^MNa dE,,J dE, ^^^E,) Pj^ {E,)a^{E,) ff{E,,,-E,) , (2 

JEi Jo L 



where Ei^^ = 0.2GeV x i. Here M is the detector mass (g), Na = 6.017 x 10^^ (mol"^) 
is the Avogadro number, is the i/a flux* — i^ij,, Pfj,, ^e-i ^e) of the T2K f^-beam [38], 
which is dominated by v^^^ but has secondary P^, v^, components. P^^^^i,^ denotes the 
neutrino oscillation probability for v^xi'^e ^ixi or P^i^, Pe ^in^e-i including the matter 
effect. ap{E^) is the cross section of the CC events for the CCQE process {X — 
CCQE) and the non-CCQE processes {X — Res) per nucleon in water. The last term of 
eq. (28), /j^(£'rec; -E'l/) is the smearing function of eq. (18) for the CCQE events, and that 
of eqs. (20) and (21) for the "resonance" events. The index D tells the detector location; 
the baseline length for D — SK is 295 km and that for the far detector D = Kr is chosen 
between L = 1000 km and 1200 km. 

The effective CCQE cross section per nucleon is slightly smaller than the naive cross 
section at high energies; 

1.0, (for £;JGeV] < 0.9) 

aJ°Q^(£;.) = ^^(£;.) x{ JK- 0.90\ ^ ^ (29) 



1-0 - 0.054 ( £ _ j , (for E,[GeV] > 0.9) 

because of occasional emission of 7r° or 7 from the oxygen nuclei. As for the naive CCQE 
cross section per nucleon, ag'"'^^{E^) for i/p {i/p = z^^, z^^, t'e, '^e) in water, we use the 
estimates of ref. [39] throughout the present analysis. The reduction factor in eq. (29) is 
our parameterization of the outputs of nuance [30]. 

The effective resonance event cross section a^^{Ei,) is the total cross section of all the 
non-CCQE CC events that satisfy the CCQE selection criteria of eq. (16). They are slightly 
different between i/^ and i/g CC events, and we find that the following parameterizations 

(0 4^^\ 
0.789 + 0.00738 log£;, - j , ( for E^GeY] > 0.51) , (30a) 

a^'^iE,) = a^^'^^iE,) (^0.810 + 0.00738 log - , ( for £;,[GeV] > 0.54) , (30b) 

reproduce well the results of nuance [30]. The gradual increase of the non-CCQE rates with 
El, reflects the growth of the number of contributing resonances and deep-inelastic events 
at high energies. 

Both the fudge factors in eqs. (29) and (30) and the smearing functions eqs. (18), (20), 
and (21), are obtained for and z/g CC events. They can be slightly different for z/^ and z/g 



tXhe flux distribution used in this report are available from the authors, or directly from the web- 
site [37]. 



CC events because of isospin breaking {nip ^ rUn, m/^+ ^ niAo, etc. ) and the presence of 
isolated protons in a water molecule. However, because the secondary anti-neutrino fluxes 
are small, we use the same fudge factors and the smearing functions for anti-neutrinos, 
simply by replacing the CCQE cross sections by those of anti-neutrinos. 
The total number of the signal CC events in each bin is now expressed as 

Ki"" = E (^m) + K5(^e) + N^iiP,) + N^iiPe)] , (31) 

X=CCQE,Res 

ior a — II and e, if there are no background. Here and £e ^-re the detection efficiencies 
for observing the //^ or signal, respectively. In actual experiments, there is a small 
probability of a percent level that a is misidentified as an e"^ signal, Pe/ij., and also the 
reciprocal probability, -P^/e, of taking as //^. In addition, significant fraction of single 7r° 
production events via NC cannot be distinguished from the CCQE signal as explained 
in the previous subsection. After adding those backgrounds the total number of a observed 
events can be expressed as 

7V;^ = (1 - Pe/^)iV;5^ + P^/e • iV^'g^ (32a) 

K,D = A/. ■ iv;5^ + (1 - p./e)iv;'g^ + iv:',rS ' (32b) 

where N^o'^ is the event numbers from the NC 7r° background in the i-th bin. 

In Fig. 5(a) and (b), typical e- and /x-like event numbers with 5 x lO^^POT for the 3.0° 
OAB at SK is shown, when the normal hierarchy is assumed. Figures 5(c) and (d) are 
for the inverted hierarchy. The histogram gives the total event numbers, and the circles 
and the triangles give the CCQE and non-CCQE "resonance" event numbers, respectively. 
The squares and the diamonds in (b) and (d) show the background event numbers from 
the misidentified vr*^ and fi"^, respectively. Events with Prec < 0.4 GeV (i = and 1) are 
not shown because we do not use them in our analysis. 

The input values of the neutrino mass and mixing parameters adopted for Fig. 5 are 

ISml^l = 2.5xlO"^eV^ sin^ ^atm = 0.5 , (33a) 
6ml2 = 8.2 X lO-^eV^ , sin^ 2^sol = 0.83 , (33b) 
sin2 2^HCT = 0.10, 5mns = 0°. (33c) 

Although the central values of the most recent measurements in eqs. (2) and (4) are slightly 
different, we use the above values in order to compare our results quantitatively with those 
of the previous studies in ref. [23, 24]. 

The matter density along the baseline between J-PARC and SK, and that between 
J-PARC and the far detector in Korea are taken as 

PSK = 2.6g/cm3 for SK , (34a) 
PKr — 3.0g/cm^ for Korea. (34b) 

These average matter densities along the baseline are obtained [34] from the recent geo- 
physical measurements [40, 41] which have typical errors of about 6%. The value for the 
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Figure 5: Typical numbers of the /i-like events, (a) and (c), and e-like events, (b) and (d), 
for the 3.0° OAB at SK with 5 x lO^^POT. (a) and (b) are for the normal hierarchy, and (c) 
and (d) are for the inverted hierarchy. The histograms gives the total event numbers, the 
circles and the triangles are the CCQE and the "resonance" event numbers, respectively. 
The squares and diamonds in (a) and (c) stand for the background event numbers from 
the misidentified 7r° and /x^, respectively. The inputs are listed in eqs. (33)-(35). We show 
only those events with Ej-ec > 0.4 GeV used in our analysis. 

T2K baseline eq. (34a) is slightly lower than 2.8g/cm^ quoted in ref. [19], because of the 
"Fossa Magna" along the baseline, in which the average density is as low as 2.5g/cm^. The 
average matter density along the baseline for the far detector in Korea depends slightly on 
the baseline length between L = 1000 km and 1200 km, because it goes through the upper 
mantle. Those details as well as the impacts of the matter profile along the baseline will 
be reported elsewhere [34]. 

Finally, the efficiencies for detecting /i^ and sinal events in eq. (31) and the probabil- 
ity of misidentifying /i^ as {Pe/fi) and that of misidentifying as /i^ (-P^/e) in eq. (32) 
are respectively. 



P. 



100% , 
1%, 



90% 



(35a) 
(35b) 



Hereafter we set P^/e = for simplicity, because P^/e ~ 1% does not affect our results 
significantly due to the smallness of the expected number of events. 



The Un survival probability is less than 40% in the region 0.4 GeV < E,^ < 1.0 GeV, 
because of the oscillation dip for Py^^^^ at ~ 0.6 GeV. Nevertheless, we expect many 
CCQE events with E^-^c < 1.0 GeV in Fig. 5(a) and (c) due to the high intensity of the 
i^fj, flux at 3.0° off-axis angle, which has a peak at E,y ^ 0.5 GeV. It catches our eyes that 
the /i-like event rate in the first bin (0.4 GeV< E^^c <0.6 GeV) is significantly larger for 
the inverted hierarchy than for the normal hierarchy. This is because the oscillation phase 
shift, the factor i?^ in eqs. (11a) and (13b), is negative for the parameters of eq. (33) so that 
the location of the dip occurs at slightly higher E^, for the inverted hierarchy. Such small 
difference in the dip location between the two hierarchies, however, can be compensated 
by a small shift in l^^riisl of several percent order. This in turn tells that |5^i3| cannot 
be measured beyond the accuracy of several percent unless the mass hierarchy pattern is 
determined; see discussions in section 5.4 for more details. 

Typical e-like events at SK are shown in Figs. 5(b) and (d). The CCQE events dominate 
the e-like events for both mass hierarchies. Because there is little high energy tail for the 
3.0° OAB and the probabihty of misidentifying 7r° as is not large at Ei, < 1.0 GeV, as can 
be seen from Figs. 4(a) and (b), respectively, the 7r° background events given by the squares 
do not dominate over the CCQE signal events. Nevertheless, they consist of about 20% 
of the total number of e-like events in the first three bins oi E^ < 1.0 GeV. Quantitative 
estimate of the 7r° background should hence be essential to measure the — > transition 
probability with confidence. 

In Fig. 6, we show the E'rec distributions of the /x-like and e-like events expected for a 
100 kton far detector at L = 1000 km and with the 0.5° OAB, for exactly the same model 
parameters of eq. (33) and the systematics of eq. (35), but with the average matter density 
of eq. (34b). 

The Ej-ec distributions of the //-like events are shown for the normal and inverted hi- 
erarchy in Figs. 6(a) and (b), respectively, where little dependence on the mass hierarchy 
pattern can be observed. The i/^ — > oscillation dip at E^, ~ 2.0 GeV is clearly seen in 
both cases, despite the contribution from the non-CCQE "resonance" events shown by the 
triangles, which has a dip at lower i?rec- 

What is most surprising in Fig. 6 is the overwhelmingly large contribution of the vr" 
background events, shown by the squares, in the e-like event distributions, both in (b) 
and (d), respectively, for the normal and the inverted hierarchies. They dominate the 
CCQE signal at low E^^c, E^ec < 1-4 GeV for the normal hierarchy and E^^c < 1-6 GeV 
for the inverted hierarchy. This is essentially because of the hard energy (broad band) 
spectrum of the 0.5° OAB, which gives rise to copious production of single events via 
the NC. Nevertheless, the CCQE event numbers supersede the 7r° background at high 
E^ec, E^ec > 1-4 GeV for the normal hierarchy, and Ej-^c, E^^c > 1-6 GeV for the inverted 
hierarchy. The significant difference in the E^cc distributions of the e-like events expected at 
a far detector, between Figs. 6(b) and (d), in contrast to the similarity of the corresponding 
distributions at SK, between Figs. 5(b) and (d), may allow us to determine the neutrino 
mass hierarchy even in the presence of the background, since the 7r° background due to 
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Figure 6: The same as Fig. 5, but for the 0.5° OAB at L = 1000 km with a 100 kton water 
Cerenkov detector. 

the NC events do not depend on the mass hierarchy. The non-CCQE "resonance" events, 
shown by the triangle, behave similarly to the CCQE signal events; the number of events 
is enhanced for the normal hierarchy and suppressed for the inverted hierarchy. Therefore, 
we expect that the contribution from the "resonance" events will enhance the sensitivity 
of the T2KK experiment to the mass hierarchy. 



4 Analysis Method 

In order to quantify the physics potential of the T2KK neutrino oscillation experiment, we 
introduce a function 

= xIk + xL + Xsys + Xjara , (36) 

which measures the sensitivity of the expected measurements on the physics parameters 
such as the neutrino mass hierarchy, sin^ 29^ct and (5mns! in the presence of statistical errors 
as well as various systematic errors including the uncertainties in the other parameters of 
the three neutrino model. 

The first two terms in eq. (36), Xsk and Xxr; respectively, measure the constraints from 



the measurements at SK and a far detector in Korea; 



E 



input 



+ 



input 



input 



input 



, {D = SK, Kr) .(37) 



Here (iV^£))™^"* ^.nd {N^ j^)"^^^^ denotes the yU-hke and e-hke event numbers, respectively, 
at SK [D = SK) and at a far detector in Korea {D = Kr), in the i-th bin of Ej-^c calculated 
as in eq. (28)- (32), and its square root gives the statistical error. The summation is over 
all bins from 0.4 GeV to 5.0 GeV at both detectors for A^^, 0.4 GeV to 1.2 GeV at SK, 
and 0.4 GeV to 2.8 GeV at Korea for N^. In order to compare our results quantitatively 
with those of the previous studies in ref. [23, 24, 25], we use the same input values of the 
neutrino model parameters, as in eq. (33), when calculating the expected number of events 
in each bin. 

The event numbers for the fit, (-/V^^o)^* and (NI^d)^^ calculated as 



(1 - Pe/.) E ^JZfpK>») 



X,a,f3 



.X,a,P 



(38a) 



, (38b) 



where the initial neutrino flavor, i/^, i/^, v^, Pg, are denoted as z/^ with a = fi, fi, e, e, 
respectively, and the superscript X denotes the event type, X = GGQE for the signal, or 
X = Res for the non-CCQE "resonance" events that pass the CCQE selection criteria of 
eq. (16). The subscript (3 distinguishes neutrinos {fi = v for or v^) and anti-neutrino 
{(3 = v for z/^ or z/g), while D = SK or D = Kr as in eq. (37). We introduce 17 normalization 
factors whose deviation from unity measures systematic uncertainties, 15 of which appear 
explicitly in eq. (38); fy for the fiducial volume and /^^ for the initial neutrino flux at = 
SK and D = Kr, for the CC cross section of X = CCQE or X = Res with neutrino 
{(3 = v) or anti-neutrino {j3 = u), and /^o for the NC cross section of producing the single 
tt" background. In addition the factor takes account of the uncertainty in the average 
matter density along the baseline between J- PARC and SK (D = SK) or Korea {D = Kr), 
which appear in the computation of the oscillation probability Pj^^^^^^ by modifying the 
matter density as 



D input 



fU „in 

Jo Pd 



(D = SK, Kr) 



(39) 



By using the above 17 normalization factors, the detection efficiencies (Sg and £^) and 



the /i-to-e misidentification probabihty (Pe/fj) , we estimate the systematic effects as follows; 



2 



Xsys 




, . , , . , , (40) 

0.05 J V 0.01 J V 0.01 ' ^ ^ 

All the errors in the first row of eq. (40) depend on the detector and its location, D — 
SK and D — Kr. The first term is the uncertainty of the fiducial volume, for which 
we assign 3% error independently for SK (fy^) and a far detector in Korea (fy^). The 
second one is for the matter density uncertainties along the T2K (Z^^) and the Tokai- 
to-Korea (f^^) basehne. The dominant source of the error in the matter density arises 
when the sound velocity data are translated into the matter density [34, 42], and we 
assign 6% error independently for each baseline. The last term of the first row is for 
the overall normalization of each neutrino fiux, which are taken independently for each 
neutrino species and the detector location. This is a conservative estimate, since it is 
likely that all the fiux normalization errors are positively correlated. The second row gives 
the uncertainty in the cross sections. Because the CCQE cross section for Ug and are 
expected to be very similar theoretically, we assign a common overall error of 3% for i/e 
and i/^ (/^°^^) and an independent 3% error for and (/^^^^). For non-CCQE 
"resonance" events (ff^), we assume 20% error ior P — u and (3 — v independently, since 
it depends not only on the single tt production cross section but also on the momentum 
distribution and the detector performance. We allow 50% error for the NC cross section of 
producing single 7r° background (/tt"), since it takes account of the uncertainty in the 7r°- 
to-e misidentification probability (P^o/g). The systematic errors in the last row of cq. (40) 
account for the performance of a water Cerenkov detector. The first and the second terms 
denote the uncertainty of the detection efficiency for e- and //-like events, respectively. In 
this analysis, we adopt = 5% and Se^ = 1%, which are taken common for SK and a 
far detector in Korea. The last one is the probability of misidentifying a /x-cvent as an 
e-event, for which a common error of 1% is assumed. In total, we adopt 20 parameters in 
simulating the systematic errors. 

Finally, Xp^ra accounts for external constraints on the model parameters: 

2 _ / {Smj^f- 8.2 X 10- / sin' 2g|^^ - 0.83^ 
^P^^ ~ 1 0.6 X 10-5 ) ^ V oW 

^ |^ sin'2e.-sin'2Cr y ^^^^ 

Although the errors of the smaller mass-squared difference and the solar mixing angle in 
eq. (41) are somewhat larger than their most recent values in eq. (4), we stick to the above 



estimates in order to compare our results quantitatively with those of the previous studies 
in ref. [23, 24, 25]. In the last term, we assume that the planned future reactor experiments 
[16, 17, 18] will measure sin^ 2^rct with the uncertainty of 0.01. 

5 Mass hierarchy 

In this section, we study the sensitivity of the T2KK experiment on the neutrino mass 
hierarchy. First, we look for the best combination of the off-axis angle at SK and the 
location of a far detector in Korea, which can be parameterized in terms of the baseline 
length L and the off axis angle from the beam center. Second, we examine carefully the 
impacts of the systematic errors, including the contribution from the uncertainty in the 
TT*^ background. In subsection 5.3, we show the sensitivity of the T2KK experiment on 
the neutrino mass hierarchy, as contour plots on the plane of sin^ 2^^rct and (5mns- In last 
subsection, we show the impacts of the mass hierarchy uncertainty on the measurement of 
\5ml^\. 

5.1 The best combination 

Here we repeat the analysis of ref. [23, 24] in which the combination of the off-axis angle 
at SK and the location of a far detector in Korea that maximizes the sensitivity to the 
neutrino mass hierarchy has been looked for, by assuming a water Cerenkov detector of 
100 kton fiducial volume at a distance between L — 1000 km and 1200 km from J-PARC. 
It should be noted here that because the detector should be placed on the earth surface, 
the allowed range of the off-axis angle at a far detector depends on the off-axis angle at SK. 
For instance, the off-axis angle observable in Korea is larger than 0.5° for the 3.0° OAB at 
SK, while it is larger than 1.0° for the 2.5° OAB at SK. 

We show in Fig. 7 the minimum Ax^ expected for the T2KK experiment after 5 x 
10^^ POT exposure as a function of the off-axis angle and the baseline length (L) of 
the far detector in Korea, when the off-axis angle is 3.0° at SK. Figure 7(a) shows the 
results when the normal hierarchy is assumed in generating the events and the inverted 
hierarchy is assumed in the fit. The opposite case, the results when the events are generated 
for the inverted hierarchy and the normal hierarchy is assumed in the fit are shown in 
Fig. 7(b). The solid-circle, open-circle, open-triangle, open-square, and open-diamond, 
give the minimum A^^ for the baseline length L = 1000 km, 1050 km, 1100 km, 1150 km, 
and 1200 km, respectively. The results depend strongly on the input values of sin^ 29j^ct 
and 6mns- sin^ 26'rct = 0.1 is assumed for all the plots and Smns is 0°, 90°, 180°, and —90°, 
from the left to the right plots. All the other input parameters are listed in eqs. (33)-(35). 
In each plot, we show by the cross symbols the highest Ax^j^ values of the previous study 
in ref. [24]. When they are higher then 30, the cross symbols are given on top of the frame 
and their values are shown in parentheses. 

All the plots in Fig. 7 confirm the trend observed in the previous studies in ref. [23, 24] 
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Figure 7: (a): Minimum Ax^ of the T2KK experiment as a function of the off-axis angle 
and the baseline length (L) of a far detector from J-PARC, after 5 x 10^^ POT exposure of 
the 3.0° OAB at SK with a water Cerenkov detector of 100 kton fiducial volume in Korea. 
The normal hierarchy is assumed in generating the events and the inverted hierarchy is 
assumed in the fit. The solid-circle, open-circle, open-triangle, open-square, and open- 
diamond, shows Ax^jjj for L = 1000 km, 1050 km, 1100 km, 1150 km, and 1200 km, 
respectively. We take sin^ 2^i,ct = 0.1 for all figures and 6mns = 0°, 90°, 180°, and -90°, 
from the left to the right plots. All the other input parameters are listed in eqs. (33)-(35). 
The cross symbols show the highest Axi^i^ value of the previous study in ref. [24]. (b): 
The same as (a), but the inverted hierarchy is assumed in generating the events and the 
normal hierarchy is assumed in the fit. 



that the sensitivity to the neutrino mass hierarchy is highest when the off-axis angle at a 
far detector is smallest and that there is little dependence on the baseline length between 
1000 km and 1200 km. This is essentially because the first oscillation maximum in the z/^- 
to-i^e transition probability occurs at around = 2 GeV in Korea, which can be observed 
via the wide-band beam of small off-axis angle but not with the narrow-band beam with 
^2.0° off-axis angle [23, 24]. It is rc-assuring that the mass hierarchy pattern can still 
be determined at Sa level just by adding a 100 kton level water Cerenkov detector at 
a right place (off-axis angle <, 1°) in Korea during the T2K experimental period (5 x 
lO^^POT), even after the realistic estimation for the reconstructed energy resolution and 
the background from single tt" production via neutral current arc taken into account. 

Unfortunately, the reduction of the Ax^^j^ values from the previous results are most 
significant at lower off-axis angles ( 1°) where the mass hierarchy discrimination power 
of the T2KK experiment is highest. This is because the high-energy tail of the wide-band 
beam that gives the high sensitivity to the mass hierarchy also gives rise to the higher rate 
of the single 7r° events via the neutral currents, as shown in Fig. 4(a). This results in the 
larger 7r° background to the i/^-to-i/g oscillation signal at a far detector; see Figs. 6(b) and 
(d). At the most favorable location of 0.5° OAB at L = 1000 km, the reduction in Axmin 
is as large as 40% to 60%, depending on 5mns and the hierarchy. We also note that the 
^MNs-dependence of the sensitivity to the mass hierarchy is somewhat smaller than that of 
the previous analysis: For instance, the reduction of A^^j^ value is largest for 5mns — 

180° 

in Fig. 7, where the highest A^^j^ value was reported in ref. [24]. This is because the 
contribution proportional to cos(5mns in the "phase-shift" term in eq. (14b) is made 
less effective in discriminating the hierarchy by the smearing in E'^ec due to the nucleon 
Fermi motion and the finite detector resolutions, which have not been taken into account 
in ref. [23, 24]. 

In Table 2, we fist the pull factors of all the parameters for systematic errors at A^^jj^, 
for 3.0° OAB at SK and 0.5° OAB at L = 1000 km, for the normal hierarchy and for all the 
four (5mns values in Fig. 7(a). It is clearly seen that the pull factors for sin^2^RCT, f^^i ^ei 
and Pe//i are most significant. The sin^ 2^^|J^^ is shifted upwards in order to compensate for 
the small event numbers expected for the inverted hierarchy. The matter density between 
J-PARC and Korea is reduced to make the matter effect in the wrong sign small. On the 
other hand, psK is slightly shifted in the positive direction, because it is the difference in the 
matter effects along the two baselines that is sensitive to the mass hierarchy. The positive 
pull factors of e^. and Pg/^t also increase the number of e-like events at a far detector in 
Korea. Reduction of these errors, in particular that of sin^ 2^^hct by the next-generation 
reactor experiments, should hence improve the sensitivity of the T2KK experiment on the 
neutrino mass hierarchy. On the other hand, the fraction of the systematic errors in the 
total Ax^-^ is not large for a 100 kton detector with 5 x lO^^POT, as shown in the bottom 
line of Table 2. Therefore, a larger detector and/or higher beam power will improve the 
sensitivity of the experiment. 

In Table 3, we show how Ax^^j^ changes from the values in ref. [24] by adding successively 





rinput 

Omns 


parameters 


0° 


90° 


180° 


-90° 


sin^ 2^RCT 


0.74 


0.18 


0.90 


1.8 


sin^ 26'soL 


0.024 


-0.010 


0.052 


0.12 




0.14 


0.067 


0.35 


0.48 


fSK 

J p 
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0.031 
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-0.13 
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fSK 
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0.24 


fKr 
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0.22 
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0.48 
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1.3 
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-^para~^-^sys 
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1.2 


4.0 


7.6 


(Xpara~l~Xsys) /^Xmin 


0.13 


0.091 


0.17 


0.27 



Table 2: The pull factors of the parameters that characterize systematic errors at A^^j^ 
for 3.0° OAB at SK and 0.5° OAB at L = 1000 km, when the normal hierarchy is assumed 
in generating the events and the inverted hierarchy is assumed in the fit. The model 
parameters are taken as in eqs. (33)-(35) for sin^ 2^r,ct = 0.1 and 5mns = 0°, 90°, 180°, and 
—90°. The pull factors whose magnitudes are larger than unity are shown by bold face 
letters. The bottom lines give the squared sum of all the pull factors and its fraction in 
the total Ax^n- 
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previous results [24] 


22.9 
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(1) 


PSK/pKr = 2.6/3.0 X (1 ± 0.06) g/cm^, = (100l°)% 


22.8 


30.3 


54.0 
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(2) 


£e = (90 ± 5)% 


20.4 


26.8 


47.4 


42.7 


(3) 


E'rec for event energy with detector resolution 


18.3 


23.3 


39.8 


37.1 


(4) 


Pe/M = (1 ± 1)% 


17.4 


19.8 


31.7 


31.5 


(5) 


7r° background 


11.1 


10.3 


20.7 


23.2 


(6) 


non-CCQE "resonance" events 


14.2 


12.7 


23.8 


28.0 



Table 3: Changes in Ax^^^ with various assumptions on the systematic effects of the T2KK 
experiment, for the combination of 3.0° OAB at SK and 0.5° OAB at L = 1000 km, when 
the normal hierarchy is assumed in generating the events and the inverted hierarchy in the 
fit. We take sva^2e^cT = 0.1 and 5mns = 0°, 90°, 180°, and -90° as inputs and all the 
others are as in eqs. (33)-(35). The top (0) row gives the previous results of ref. [24], and 
each row gives the results after changing the conditions one by one. The bottom (6) row 
gives our results shown in Fig. 7(a) 

the effects introduced in this analysis, for the combination of 3.0° OAB at SK and 0.5° 
OAB at L = 1000 km, when the normal hierarchy is assumed in generating the events 
and the inverted hierarchy is assumed in the fit. The first row (0) gives the results of the 
previous study in ref. [24]. In the row (1), we change the average matter density along the 
T2K baseline from 2.8 to 2.6 g/cm^ and the error of psK and pxr are doubled from 3% 
to 6%, and we also introduced a 1% error in the ji detection efficiency. The Ax^ values 
are slightly reduced for 5mns = 180° and —90° cases, mainly because of the increase in the 
matter density errors. In the row (2), we further introduce the detection efficiency for the 
e-like events, Se = (90 ± 5)%, and the Ax^in fo^ all 5mns decrease by about 10% reflecting 
the 10% decrease of the signal events. In the row (3), we introduce smearing in E^cc due 
to the nuclear Fcrmion motion and realistic energy resolution of detectors. Because the 
matter effects in the phase-shift term is diluted by the smearing, the decrease in AXmin 
is largest 6 = 180°; see the term proportional to cos^^ns in cq. (14b). In the row (4), we 
take into account the particle misidentification probability P^/^ = (1 ± 1)%. Since this 
change makes the fake e-like events around the dip of the z/g transition probability, 

the reduction in Axmin significant even for 1% misidentification probability, if its error 
is as large as 100%. In the row (5), the single 7r° events reduce the physics potential of the 
T2KK experiment significantly, because the — > i/g signal at small E'rec is dominated by 
the 7r° background at a far detector in Korea, as shown in Fig. 6. In the bottom row (6), 
we add the non-CCQE "resonance" events in the analysis. These events make Ax^in l^^ge, 
because their magnitudes are also proportional to the — > transition probability. 

In Fig. 8, we show the minimum Ax^, the mass hierarchy discrimination power of the 
T2KK experiment, when the beam center is 2.5° below the SK. All the other contents of 



25 



16 



! Smns=Q° : : 

1000km 
1 050km 
1100km 
□ 1150km 
O 1200km 
X Previous 



O 
A 



X 



X 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



(a) normal hierarchy (OAB:2.5@SK) 

X (48.5) 











5mns=v 


0° 












































X 




















































































X 
























• 






































■ X 



















X 



5mns- 



X 



180' 



§ ^ X 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



(43.3) 



X 



90° 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



25 



16 



<3 



X 



: '. : 3MiMs=Q° : : 

1000km 

O 1050km 
A 1100km 

□ 1 1 50km 

X Previous 



X 



m 



X 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



(b) inverted hierarchy (OAB:2.5@SK) 

X (43.8) 



X 



X 



=90° 



X 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



X 



i 



X 



3mns- 



180' 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



X 



Sun 



X 



^ B 



90°: 



0.5 1.0 1.5 2.0 2.5 3.0 
off-axis angle at Kr. (degree) 



Figure 8: The same as Figure 7, but with 2.5° OAB at SK. 



Fig. 8 are the same as those of Fig. 7. Because of the geological constraint, the 2.5° OAB 
at SK cannot provide 0.5° OAB in Korean peninsula [23, 24]. When the off-axis angle is 
2.5° at SK, the optimum OAB for a far detector in Korea is 1.0° at L = 1000 km. The 
value of '^Xmin is '^o^ significantly different between the 3.0° OAB at SK and the 2.5° OAB 
at SK, when the off-axis angle in Korea is fixed as 1.0°. It confirms our understanding 
that the energy profile or the hardness of the neutrino beam observed at a far detector is 
essential for the mass hierarchy discrimination. 

5.2 Uncertainty of the tt^ background 

In this subsection, we examine the impacts of the vr" background in more detail. In our 
analysis, we adopt the following uncertainties for the relevant cross sections 

,CCQE^^^Q^Q3^ /f^- = l±0.2, ^0 = 1 ±0.5, (42) 
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Figure 9: Minimum A^^ of the T2KK experiment as a function of the minimum -Erec, when 
only those e-hke events with larger i^rec values are retained in the analysis at a far detector. 
The events are calculated for the 3.0° OAB at SK and a 100 kton water Cerenkov detector 
at 0.5° OAB and L = 1000 km with 5 x 10^^ POT exposure, for sin^ 2^i,ct = 0.1 and 
^MNs = 0° and the parameters of eqs. (33)- (35), for the normal hierarchy, while the inverted 
hierarchy is assumed in the fit. The dashed-dotted line with diamonds, the dotted line 
with up-triangles, the solid line with circles, and the dashed line with down-triangles give 
^Xmin when the uncertainty in the 7r° background rate A/^o is 10%, 30%, 50%, and 70%, 
respectively. 

where (3 = v and z/; see eq. (40). The 3% error in the CCQE cross sections should be 
achieved in the near future, whereas there is a possibility that the non-CCQE "resonance" 
cross sections and the neutral current single 7r° production cross section can be measured 
more accurately than 20% and 50%, respectively, assumed in this analysis. We therefore 
repeat the fit by varying Aff"^ between 10% and 30%, and A/^o between 10% and 70%. 
We find little impacts of those variations on the magnitude of Ax^j^, which conform with 
the small pull factors for these parameters in Table 2. It turns out that the uncertainty 
in the non-CCQE cross section does not affect the mass hierarchy sensitivity of the T2KK 
experiment because it tends to cancel in the ratio of the //-like and e-like events. 

In case of the 7r° background to the e-like events, however, the smallness of the impacts 
of varying A/ttO between 10% and 70% is striking, and we examine the cause carefully. 
In Fig. 9, we show Ax^j^ of the T2KK experiment as a function of the lowest E^^c above 
which the e-like events are counted at the far detector in Korea. All the other conditions 
and the input parameters are the same as those of Fig. 7(a) and Table 3, for 5mns = 0°. 
The dash-dotted line with diamonds, the dotted line with upper triangles, the solid line 
with circles, and the dashed line with lower triangles are obtained with the 7r° background 
normalization error of A/^o = 10%, 30%, 50%, and 70%, respectively. 

It is clearly seen that there is little dependence on the error A/^o when we use all the 
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Figure 10: The Ax^j^ contour plot for the capabihty of the T2KK experiment to determine 



the neutrino mass hierarchy on the sin and 6mns plain, (al) and (bl) are for the 

normal hierarchy, while (a2) and (b2) are for the inverted hierarchy. The OAB combination 
for (al) and (a2) is 3.0° OAB at SK and 0.5° OAB at L = 1000 km, 2.5° OAB at SK and 
1.0° OAB at L = 1000 km is for (bl) and (b2). All the input parameters other than 
sin^ 26'"cT* and 5mns* are the same as those in Figs. 7 and 8. 

data with i^rec > 0.4 GeV as has been assumed in our analysis. As the i^rec threshold is 
increased, however, the reduction in Ax^j^ becomes significant as A/^o increases. This is 
because the normalization of the vr" background can be determined by the e-like event rate 
at low Ej-ec where the vr^ background dominates the oscillation signal; see Fig. 6(b) and 
(d). This suggests strongly that we should understand not only the overall normalization 
of the 7r° background but also the energy and angular distribution of singly produced vr^'s 
in the neutral current events as well as the momentum dependence of the error of the 
7r°-to-e misidentification probability Ps/ttO, whose parameterization is given in Fig. 4(b). 
Detailed studies of the normalization and the shape of the background should be the 
most important task before the physics case of the T2KK experiment can be established. 



5.3 Dependence of the OAB at SK 

In Figs. 7 and 8, we find that the best location of the far detector to determine the neutrino 
mass hierarchy is at L = 1000 km away from the J-PARC, where 0.5° OAB can be observed 
for the 3.0° OAB at SK (Fig. 7), or 1.0° OAB for the 2.5° OAB at SK (Fig. 8). In this 
subsection, we compare carefully the two combinations since they can be interchanged, or 
the (1.0 - 9)° OAB can be observed for the (2.5 + 9)° OAB at SK, simply by adjusting 
the beam direction at J-PARC (up to \9\ ^0.5°) for a fixed far detector location along the 
baseline at L ^ 1000 km. 



In Fig. 10, we show the contours for /^Xmin — 4, 9, 16, 25 in the plane of sin^ 2^rct* and 
5mns*- The wrong hierarchy can be excluded with the n-cr confidence level, if the true values 
of sin^ 20j^cT and 5mns lie in the right-hand side of the A^^jjj = contour. The upper 
figures (al) and (a2) are for 3.0° OAB at SK with 0.5° OAB at L = 1000 km, and the lower 
figures (bl) and (b2) are for the 2.5° OAB at SK with 1.0° OAB also at L = 1000 km. 

It is clearly seen from the figures that the mass hierarchy can be determined better by 
the combination of 3.0° OAB at SK and 0.5° OAB at L = 1000 km than the combination of 
2.5° and 1.0° for all the input values of sin^ 26'j^ct and 5mns and for both hierarchy patterns. 
For instance, by comparing the figures (al) and (bl) we find that the normal hierarchy can 
be established at 3(T level, Ax^j^ > 9, when sin^ 2^rct ^0.08 (0.12) for the combination 
of 3.0° and 0.5° (2.5° and 1.0°). Likewise, from the figures (a2) and (b2), the inverted 
hierarchy can be established when sin^2^RCT ^0.09 (0.12) for the combination of 3.0° and 
0.5° (2.5° and 1.0°). The difference is significant when |5mns| ^90° where it is difficult to 
determine the mass hierarchy. On the other hand, we find little dependence on the off-axis 
angle between 3.0° and 2.5° at SK when 5mns — 180°, where the mass hierarchy can be 
determined with relative ease. 

The reason for the strong dependence on the off-axis angle when |5mns| ^90° can be 
explained by the hardness of the 0.5° OAB that provides sufficient flux at the z/^ — > oscil- 
lation maximum around E^, ~ 2.0 GeV. It is essentially the mass hierarchy dependence of 
the amplitude shift term, in eqs. (12) and (14a), which contribute to the determination, 
and the hardness of the 0.5° OAB helps enhancing the signal. When 5mns — 180°, in addi- 
tion to the amplitude shift term, the phase-shift term S^, in eq. (12), becomes significant 
because the leading term and the sub-leading term in eq. (14b) adds up to make \B^\ large 
at cos(5mns — —1. The mass hierarchy dependence due to the phase shift term B'^ turns 
out to give significant difference in the — > transition probability at lower E^, [23, 24] , 
and the downward shift of the flux maximum E^, in the 1.0° OAB can be compensated for. 

In the absence of a concrete evidence that the nature chooses cos 5mns — 180°, it is clear 
that the effort to make the off-axis angle at the far detector as small as possible should be 
valuable. The sensitivity difference between 0.5° OAB and 1.0° OAB in Fig. 10 corresponds 
to about a factor of two difference in the product of the flducial volume of the far detector 
and the POT, the beam power times the running period. 

5.4 Impacts on the [(^m^gl measurement 

In this subsection, we comment on the implication of the mass hierarchy uncertainty in 
the measurement of the absolute value of the larger mass-squared difference. 

In Fig. 11, we show the minimum A^^ of the T2KK experiment as a function of l^m^gl 
with the optimum OAB combination of 3.0° at SK and 0.5° at L = 1000 km. Fig. 11(a) 
is for the normal hierarchy and Fig. 11(b) is for the inverted hierarchy. The flve curves 
are for sm^2e'^ST = 0.02, 0.04, 0.06, 0.08, and 0.1, which are denoted by the solid, long- 
dashed, short-dashed, dotted, and the sohd line again, respectively. The CP phase is flxed 
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Figure 11: Ax^j^^ of the T2KK experiment as a function of l^mjgl for 3.0° OAB at SK 
and 0.5° OAB at L = 1000 km. (a): The normal hierarchy case at 5mns* = 0°. The 
sohd hne, long-dashed, short-dashed, dotted, and the sohd line again, show the results for 
sin^26'™(?T* = 0.02, 0.04, 0.06, 0.08, and 0.1, respectively. All the other input parameters 
are those in eqs. (33)-(35). The blue lines, which are almost degenerate, are obtained when 
the right hierarchy is chosen in the fit, whereas the red lines are obtained with the wrong 
hierarchy, (b): The inverted hierarchy case. 



at 5mns = 0° and all the other parameters are those of eqs. (33)-(35). In both cases there 
is a set of five curves with Ax^j^ = 0° at I^J^i^gl = 2.5 x 10~^ eV^, the input value. All 
the five curves are almost degenerate in the set, which exhibits the insensitivity of the 
t'/i ~^ ^tJ- survival probability on sin lOj^Qr^'i see eqs. (11a) and (13). On the other hand, 
there is another set of five curves with Ax^j^ at 0.1 x 10^^ eV^ smaller (larger) than the 
input value when the mass hierarchy is normal (inverted) . These curves with sin^ 2^j^ct 
dependent Ax^j^^ are obtained when the opposite hierarchy is assumed in the fit. 

The larger mass-squared difference is determined from the T2KK experiment correctly 

as 

(2.5 ± 0.02) X 10"W , (43) 

if we know the mass hierarchy pattern. However, if we do not know the mass hierarchy 
pattern, the other solution 



\5mU 



(2.4 ± 0.02) X 10 ^eV^ (for the normal hierarchy) , 
(2.6 ± 0.02) X 10~^eV^ (for the inverted hierarchy) 



(44) 



appears for every sin^ 2^"^*- The wrong solution (44) are about 3.5cr away from the correct 
solution (43). The difference of =f0.1 x 10~^ eV^ in the mean value can be explained by the 
phase shift term in the —>■ survival probability; see eqs. (11a) and (13b). From 
the peak location at 

+ B' = \, (45) 



Ai3 



the location of the solution with the wrong hierarchy can be estimated as 





N fit 






\ input 




) 






) 



fo.U- 0.023 ('-^Y\osS-^\ X 10-3eV^ 

(5^2^)input \ \ 0.10 ' 



(46) 



The magnitude of the difference is almost 0.09 x lO'^ eV^ for sin^ 2^1^^"* = 0.1 at cos S'^ST = 
1, and it grows to 0.11 x 10"^ eV^ as sin^ 2^h,ct decreases, as can be observed from the 
figures. This result suggests that the absolute value of the larger mass-squared difference 
cannot be determined uniquely, if the mass hierarchy pattern is not known. Because the 
T2KK experiment can determine the mass hierarchy from the z/^ — >■ i/e transition rates 
for sufficiently large sin^26'RCT*; the fake If^m^sl can be excluded for larger sm^29^^c? as 
shown by the Ax^jj^ values of the wrong solutions in Figs. 11(a) and (b), which grow with 
increasing sin^ 2^r°(?t*. If we do not make use of the z/^ — > z/g transition signal in the fit, 
all the solutions with the wrong mass hierarchy has A^^j^ ~ 0, indistinguishable from the 
correct solution. 

Let us note in passing that the T2K experiment suffers from the same uncertainty in the 
measurement of |(5m^3|. If we drop all the data from the far detector in the above analysis 
with Sml^ = ±2.5 x 10"^ eV^, we find the fake solution with 

\5ml^\ ~ (2.4 ± 0.04) x lO^^eV^ (for the normal hierarchy) , 

~ (2.6 ± 0.04) X lO^^eV^ (for the inverted hierarchy) , (47) 

instead of eq. (44). The ^Xmin values for the wrong solutions are indistinguishable from 
zero for all the sin^ ^Oj^qt input values. The difference of about 0.1 x 10"^ eV^ between the 
correct and the wrong solutions remains the same, because the formulae (45) and (46) arc 
valid near the oscillation maximum at all baseline length L as long as the earth matter 
effect remains a small perturbation as in eqs. (11) and (13). Since the two solutions are 
about 2(7 away, the experiment should present two values of l^m^gl until the mass hierarchy 
is determined. 



6 CP phase 

In this section, we study the capability of the T2KK experiment for measuring the leptonic 
CP phase Smns with the optimum OAB combination, 3.0° OAB at SK and 0.5° OAB at 
L = 1000 km, with 5 x lO^^ POT exposure. 

In Fig. 12, we show Ax^ contours in the plane of sin^2^^jiCT and 5mns when the mass 
hierarchy is normal (mg — ml > 0). The input values arc sin 29'^^:^' = 0.10 and 0.04, and 
^MNs* = 0°, 90°, 180°, and —90° in Figs, (a), (b), (c), and (d), respectively. The other input 
parameters are those in eqs. (33)- (35). The contours for Ax^ = 1, 4, and 9 are shown by 
the solid, dashed, and dotted lines, respectively. The thick red lines show the Ax^ contours 
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Figure 12: The Ax^ contour plot for the T2KK experiment in the plane of sin^ ^O^ct and 
^MNs when the mass hierarchy is normal. Allowed regions in the plane of sin^ 29^ct and 
6mns are shown for the combination of 3.0° OAB at SK and 0.5° OAB at L = 1000 km with 
5 X 10^1 POT. The input value of sin^ 2^h,ct is 0.10 and 0.04 for 5mns = 0° (a), 6mns = 90° 
(b), 5mns = 180° (c), and 5mns = —90° (d), and the other input parameters are listed 
in eqs. (33)-(35). The input points are indicated as the solid blobs. The contours for 
Ax^ = 1, 4, and 9 are shown by the solid, dashed, and dotted lines, respectively. The 
thick red contours are obtained when the right hierarchy is assumed in the fit, whereas the 
thin blue contours with the local minimum by the solid square show the results when the 
opposite mass hierarchy is assumed in the fit. 
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Figure 13: The same as Fig. 12, but when the mass hierarchy is inverted (mg — ml < 0). 



when the right hierarchy is chosen in the fit, whereas the thin blue hues show the results 
when the wrong hierarchy is assumed in the fit. The solid blobs in each figure denote the 
input points [Ax^ = 0) and the solid squares are the local minima for the fit with the 
wrong hierarchy. 

We find from the figures that 6mns can be constrained to about ±45° at la level for 
all the input value of 5mnT and sin^2^HCT- The insensitivity of the measurement error 
of 5mns on sin^ [23, 24, 29] persists. However, at 3a level, the contour closes only 

for 5mns* = ±90° at sin^ 26^^"* = 0.1; (b) and (d). Moreover, there appears a shadow 
island where the inverted hierarchy is assumed in the fit, for all the four 6mnT cases at 
sin^ 26^"^ * = 0.04. The shadow contours cover the whole 5mns region when 5mnT = 90°, 
where the number of the z/g signal events is the smallest among the four ^mns* cases. 

These observations are in sharp contrast with the previous ones, shown in e.g., Fig. 8 
of ref . [24] , where it has been shown that the 5mns can be constrained to about ±30° for all 
the four input 6mns values at sin^ 26^"^* ^ 0.02 and that the shadow islands from the wrong 
hierarchy solution are small and they appear only for 5mns* = 90° at sin^ 26*^"* = 0.04 and 
for = 0°, 90°, and 180° at sin^ 2Cct* = 0.02. We find that both the reduction of the 
sensitivity from ±30° to ±45° and the appearance of the big shadow islands are mainly 
due to the vr" background for the e-like events, while the smearing effects due to nuclear 
Fermi motion and the detector resolution also contribute at the sub-leading level. 



In Fig. 13, we show the same contour plots as in Fig. 12, but for the inverted hierarchy 
case. We find that the la constraints on S^ns are shghtly worse than those of the normal 
hierarchy case in Fig. 12: The la error remains at about ±45° for = 0° and 180°, but 
it grows to about ±60° or larger for 6mnT = ±90°. As in the case of the normal hierarchy, 
^MNs can be constrained at 3a level only for (5mns* = ±90° at sin^ 2fe"* = 0.1. The 2cT-level 
shadow islands appear for all the four SmnT cases at sin^ 29'^^^' = 0.04, which is consistent 
with the observation of Fig. 10 (a2) where all 5mns* points lie below the A^^^^j^ = 4 contours. 
The 3(7 contours of the wrong solutions, denoted by the thin blue dotted lines, cover the 
whole 6mns region for = 0° and -90° at sin^ 2^^,^?^* = 0.04. The significant loss of 

the sensitivity to (5mns as compared to Fig. 9 of ref. [24] can also be explained by the tt" 
background to the — >■ z/g oscillation signal. 

In summary, the capability of the T2KK experiment to measure the CP phase of the 
lepton fiavor mixing (MNS) matrix is significantly worsened by the 7r° background in both 
normal and inverted hierarchy cases. This is because the large 7r° background to the 
J^n ~^ oscillation signal at the far detector, as shown in Figs. 6(b) and (d), reduce 
significantly the sensitivity to the amplitude-shift term and the phase-shift term 
which have contributions proportional to sin^MNs ^-nd cos^mns; respectively [23, 24]. These 
terms proportional to A12 in eq. (14) can be measured by comparing the shifts at a near 
(L ~ 300 km) and a far (L ~ 1000 km) detectors [23, 24, 29] without using the beam. 
Since the 7r° background worsens the measurements of and at the far detector, the 
sensitivity to Smns deteriorates significantly. The use of beam in addition to the i/^ beam 
[28, 29] may be helpful in recovering the sensitivity, since at least the detector-dependent 
errors of the tt'^ background events should be common for both beams. 



7 Summary and conclusion 

In this paper, we elaborate the previous analyses of ref. [23, 24] on the physics potential 
of the T2KK experiment by taking into account the smearing of reconstructed neutrino 
energy due to the Fermi motion of the target nucleus and the finite resolution of and 
momenta in a water Cerenkov detector. We also include the events from the non- 
CCQE "resonance" events that survive the CCQE event selection cut of eq. (16), and 
the contribution from the single 7r° production via the neutral current interactions, which 
mimic the appearance signal in a water Cerenkov detector. 

In order to estimate the reconstructed energy (-Ercc) distribution efficiently, we introduce 
the smearing functions for the CCQE and non-CCQE "resonance" events that map the 
incoming neutrino energy E^, onto the reconstriictcd energy Ej-ec by using the Mote Carlo 
event generator nuance [30]. The effect of the detector resolution for and /i^, see Table 1, 
has also been taken into account. The smearing functions for the CCQE events are given 
in eq. (18) with cqs. (A7)-(A6) for u^, and cqs. (A10)-(A9) for z/g . Those for non-CCQE 
"resonance" events are parameterized as in eq. (20) with eqs. (A15)-(A20) in the region of 



0.55 < < 1.2 GeV and eq. (21) with eqs. (A22)-(A27) for 1.2 GeV< < 6.0 GeV. 
For estimating the background from the single 7r° production, we generate single tt^ events 
from the NC interactions for each off-axis beam (OAB) also by using nuance [30], and 
parameterize the probability that a 7r° is misidentified as an e^-like event, Pe/w, in terms 
of the energy ratio and the opening angle of the two photons for the 7r° decay- in-flight; see 
Fig. 4(b) and eqs. (26) and (27). 

We study the sensitivity of the T2KK experiment on the neutrino mass hierarchy by 
placing a water Cerenkov detector with 100 kton fiducial volume at various location in 
Korea for the 3.0° and 2.5° OAB at SK. The neutrino beam at an off-axis angle greater than 
about 0.5 (1.0) can be observed in Korea, at the baseline length 1000 km <^L <^ 1200 km, 
for the 3.0° OAB (2.5° OAB) at SK. We find that the highest sensitivity is achieved 
for the combination of 3.0° OAB at SK and 0.5° OAB at L = 1000 km, confirming the 
results of ref. [23, 24]. With 5 x lO^^ POT, which is the planned exposure of the T2K 
experiment, the mass hierarchy can be determined at 3(T level if sin^ 2^jj^ct ^ 0.08 (0.09) for 
the above OAB combination, when the neutrino mass hierarchy is normal (inverted). For 
the combination of 2.5° OAB at SK and 1.0° OAB at L = 1000 km, the 3(7 sensitivity is 
obtained for sin^ 26'rct ^0.12 for both hierarchies; see Fig. 10 in section 5.3. These figures 
show significant reduction of the sensitivity as compared to the results of the previous 
studies, such as Fig. 6 of ref. [24], which show that the neutrino mass hierarchy can be 
determined for sin^ '2.9^ct ^0.05(0.06) at 3(7, when the hierarchy is normal (inverted), with 
the same combinations of the OAB's, and with the same detector size and the POT. 

We find that the main cause of the reduction in the sensitivity is the background from 
the single 7r° production; see Table 3 in section 5.1. The smearing in the reconstructed 
energy has a significant effect when (5mns — 180°, where the mass hierarchy dependent 
oscillation phase-shift term is large. The contribution from the non-CCQE "resonance" 
events help discriminating the mass hierarchy, because these events are also a part of the 

~^ oscillation signal. 

We also examine the prospect of the CP phase measurement for the T2KK experiment 
with the above OAB combination. The sensitivity of the 5mns measurement is also reduced 
significantly from that of the previous study in ref. [24], which found the la error of about 
±30°, to about ±45° or even ±60° in some cases. The main cause of the worsening of 
the error is again the tt" background for the e-like events at the far detector that makes it 
difficult to measure the baseline dependence of the — >■ Vf. oscillation amplitude and the 
phase: sIu^mns is measured by the amplitude difference and cos(Jmns is measured by the 
phase difference [23, 24]. 

The 7r° background reduces significantly the physics potential for the mass hierarchy 
determination and the CP phase measurement of the T2KK experiment. If we understand 
better the physics of the -k^ production and its decay signal inside the water Cerenkov detec- 
tor, the sensitivity of the experiment on these fundamental parameters should be improved. 
Detailed investigation of the normalization and the shape of the 7r° background should be 
one of the most important tasks to evaluate quantitatively the physics discovery potential 



of the T2KK experiment. 
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Appendix A Smearing functions /^{Ej-qciEi,) 

In the appendix, we show our parameterization of the smearing functions, f^[E^f.c'-,Ev), 
which map the incoming neutrino energy, E^^ onto the reconstructed energy, E^cci for the 
quasi-clastic events. The superscript X denotes the event type, X = CCQE for the CCQE 
events, or X = Res for the non-CCQE "resonance" events that pass the CCQE selection 
criteria of eq. (16), and the subscript a is for /i or e: a = fi for z/^ events and a = e for Ue 
events. These functions take account of the Fermi motion of the target nuclcon inside the 
oxygen nucleus and the finite energy-momentum resolutions of a muon and an electron in 
a water Cerenkov detector listed in Table 1. 

Appendix A.l CCQE events 

The Erec distribution of the CCQE events, which are generated by nuance [30], can be 
parameterized accurately by 3 Gaussians, 

L (ErecE^) -^;^^}^JniE.)ew 2«(i?.))2 ) ' 

in the region of 

0.3 GeV < < 6.0 GeV for 0.4 GeV < £;rec < 5.0 GeV . (A2) 
The index a takes /i for i/^ and e for i/g events, and each function is normalized by 

3 

A-{E,) = V^^<(£;.)<(£;.) . (A3) 

n=l 



The variance a", the energy shift SE". and the normahzation factors rfa are functions of 
the incoming neutrino energy E,^, with r"{Ey) = 1. The first and the second Gaussians 
account mainly for the nuclear Fermi motion, and we can set 6E2 = 5E^. The third 
Gaussian is necessary to account for the asymmetry in the E,.^ — E^ distribution such 
as the Fermi block effect at the low energies, and the asymmetric momentum resolution 
effects at high energies. 

All the coefficients are parameterized compactly by using the variables 



X = [GeV] - 1 , ^ = /MGeV]-l, 

which vanish at = 1 GeV. The variance cr^(MeV) of the three Gaussians 

= 39.7 + 68.5^, 



(To 



cx5 



82.7 - 50.2a; + 259^, 

197 + 486a; - 606^ + 203a;C , 



(A4) 



(A5a) 
(A5b) 
(A5c) 



and the energy shift terms 5E!^ (MeV) 



5Ef 



5Ei^ = 
SE^ = 



35-2.5a;(l-1.2a;)/(l + a;), 

■l + 5.27x- 8.67a;2 



0.053 + 0.033X 



1.83a;=^ 



1 + 3.65x + 4.35a;2 
are given in units of MeV. The normalization factors are — 1 and 

= 1.1- 0.96a; + 0.44a;^- 0.076a;^ + 0.0047a;^ 



(A6a) 
, (A6b) 

(A7a) 

= 0.365 - 1.97^(1 - 0.634^ + 0.464^^ + 0.293^^ - 0.342^^)/(l + 1.23a;) . (A7b) 

The first and second variances are determined mainly by the sum of the nuclear Fermi 
motion and the momentum resolution of the water Cerenkov detector. In the absence of the 
momentum resolution error, two Gaussians, one with a constant variance of ~ 60 MeV and 
the other with a larger variance of ~ 190 MeV at ~ 1 GeV which decreases slowly with 
energy, can account for the bulk of the Fermi motion effects on the E^ec — E^, distribution; 
see Fig. 1. It is the smearing effect due to the energy resolution which increases the first 
two variances as \/E^ at high energies. The value of 2 does not depend on E^ much, 
because they are essentially determined by the nucleon and lepton masses; see eq. (15). The 
third Gaussian has much larger variance than the first two, and it accounts for the Fermi- 
blocking effect at small E^ and the momentum resolution asymmetry at high energies. 
Consequently, rg is significant only at low energies [E^, < 0.7 GeV) and at high energies 
{E^ > 4 GeV). 

For the u^, case, the variance cr^(MeV) is expressed as 



^1 



55.5 - 19.6a; + 98.9^ , 

125- 51.3a; + 201^, 

273 - 102a; + 1560{ + llla;^. 



(A8a) 

(A8b) 
(A8c) 



the shift term SE^{MeY) is given as 

5EI = 5E^ = 40 - 0.99x + 3.3x^ - O.llx^ - 2.2x/(l + x) , 

SEI = (T^ [-0.16 + 0.68a; - 2.6C + l.lx/(l + x)] , (A9a) 

and the normahzed factors are rf — 1 and 

= 0.67 - 0.58x + 0.58x^ - 0.16x^ + 0.019x\ (AlOa) 

r| = 0.094 - 0.040a; + 0.031x2 - 0.016a;^ + O.OOSQx^ (AlOb) 

The three variances in eq. (A8) behave similarly to those for z/^, but cr^ is larger than 
o"^, because the energy resolution of the e-like events are worse than that of the /i-like 
events; see Table 1 in section 3. The energy shifts SEI2 behave similarly to 5Ei2) while 
5 El differs significantly from ^£'3 at low energies, because the asymmetry of the E-^^c — 
distribution in the sub-GeV region is sensitive to the mass and the momentum resolution 
of the emitted charged lepton. The normalizations g behave similarly to rg g, except at 
very low energies {E^, ^0.7 GeV) when the muon mass in not negligible and at very high 
energies {E^, ^ 3.5 GeV) due to resolution effects. 

Appendix A. 2 Nuclear Resonance Events 

The Ey-cc distribution generated by nuance [30] for the non-CCQE events that pass the 
CCQE selection criteria of eq. (16) is also parameterized for i/^ or Ug. We find that 3 
Gaussians 

fnE,.,E. < 1.2 GeV) = ^t.„"(^„,exp , (All) 

suffice in the region of 

0.55 GeV <E^< 1.2 GeV for 0.4 GeV < ^^ec < 5.0 GeV , (A12) 
whereas 4 Gaussians 

/r(^.; E.>V2 GeV) ^ ^ t f:(E, exp (- . (A13) 

are necessary in the region of 

1.2 GeV <E^< 6.0 GeV for 0.4 GeV < E,^^ < 5.0 GeV, (A14) 

because the number of contributing resonances grow at high energies. Here again a — /i 
for i/^ and a = e for i/g events, and the functions are normalized as in eq. (A3). 

By using the same variables x and ^ in eq. (A4), the variances (T^(MeV) of the 3 
Gaussians in eq. (All) are 

o-f = 97.8 - 3670a; + 7410^ + 1540a;^ , (Al5a) 
= 98 + 390a; + 500a;^ (Al5b) 

= 27 - 22a; - 33a;2 - 480a;=^ - 1000a;^ (A15c) 



and the energy shift terms 5£'^(MeV) are 



= 382 -8170a; + 16200^ + 3670a;^, (A16a) 
SE^ = 579 - 22600X + 45800^ + 10800x^ , (Al6b) 
5E^ = 210 - 12x - 7.8a;2 - 170a;^ (Al6c) 

and the normalized factors are = 1 and 

= 0.42 + 0.16X- 1.9x^ (Al7a) 
= 0.1 - 0.202a;/ (1 + 2.03a;) . (Al7b) 

The first Gaussian is mainly related to the A-resonance. The order of the first variance 
is similar to the sum of the width of the A(1232), Fa — 60 MeV [15], the Fermi motion of 
the target, oTcrmi — 60 MeV, and the momentum resolution asp/p ~ 30 MeV at E',^ ~ 1 GeV. 
The value of SE^'^ is roughly the distance between the peak of the CCQE events and that 
of the A events, which is about 400 MeV. The second Gaussian with growing variance 
of about 100 MeV at -Ej, = 1 GeV and with larger energy shift of 600 MeV accounts for 
contribution of N(1440) and higher resonances. The third Gaussian is necessary to take 
account of the nuclear effects and the asymmetry from the momentum resolution. 

For the Ue case, we find 



and 



SE: 

both in MeV units, f ^ = 1 and 



al = 102 - 2960a; + 5940^ + ISOOx^ , (Al8a) 
= 114 + 370.0a; + 244x^ (A18b) 
= 27.1 - 717a; + 1390C + 425a;^, (Al8c) 



SEt = 372 - 9530x + 18900^ + 4630a;^ , (Al9a) 
5E^ = 580.0 - 15600a; + 31400^ + 6070a;C , (Al9b) 
= al (7.6 + 5.3a; - 9Ax'^ - 24a;^) , (Al9c) 



= 0.37 - 5.8a; + 13^ , (A20a) 
fl = 0.067 - 0.068a; - 0.37a;^ - 3.4a;^ (A20b) 



There are no large difference between Ue and for all variances, energy shift terms, 
and the normalization factors. The small differences are mainly due to the difference in 
the e and /i momentum resolutions. 

In the high-energy region of eq. (A 14), we introduce variables y and rj 



2/ = [GeV] - 2 , V = VEu [GeV]/2-l, (A21) 



which vanish ect E,^ — 2 GeV. For i/^ {a — fx), the four variances are 

5-f = 110 + 12^/ + ^/^ (A22a) 

= 160 + 36?/ + 2. V, (A22b) 

5-^ = 320 - 1240?/ + 5640r7 + 328|/77 , (A22c) 

= 177 + 2580y - 91 907^ - 747^/77 , (A22d) 

and the energy shift terms are 

5E^ = 350 - 1.2y + 1.4y^ (A23a) 

SE^ = 530 + 20y - 4.9y2 + 0.58y^ (A23b) 

= 823 - (3550y - I4IOO77 - ISAOyrj) / (1 + 0.921y) , (A23c) 

SE^ = 1500 + 990y - 1500?7 - 760|/?] , (A23d) 

both in MeV units. The normalization factors are — 1 and 

= 0.31-0.18y + 0.22y2- 0.096y^ + 0.0152/\ (A24a) 

= 0.194 + 0.316?/- 1.3777-0.1221/77, (A24b) 

= 0.0369 - 0.07287/ + 0.25477 + 0.03567/77 . (A24c) 

The first Gaussian is mainly related to the A-resonancc; ~ 100 MeV and 6E!^ ~ 



350 MeV at all energies, similarly to a'^ and 5E'^ in eqs. (A15) and (A16). Because the 
number of the resonance modes which contribute to the second Gaussian increases with Ei,, 
the second variance grows from ~ 160 MeV at = 2 GeV to 240 MeV at = 4 GeV. 
A'"(1440) dominates the second Gaussian, and 5E^ does not grow much from 530MeV at 
E^ — 2 GeV. The resonances of mass greater than 2 GeV contribute to the third Gaussians; 

~ 320 MeV and SE^ ~ 830 MeV at = 2 GeV, which grow to 450 MeV and 900 MeV, 
respectively, at £'1, = 4 GeV. The last Gaussian is necessary to reproduce the tail at low 
energies. 

For the Ue events at £'1, > 1.2 GeV, we find 

at{y) = 110 + 147/ - 0.507/2, 

a^{y) = 181 + 2 1 507/ + 2077/2 - 844077 - 18907/77, 

al{y) = 334 - 17507/ + 762077 + 5567/77, 

al{y) = 222 + 103007/ - 403077 - 3570y77 , 

and 



5El{y) = 360 - y - 1.6y2 + 0.95y^ - 0.026yS (A26a) 

5E^{y) = 510-61.1y + 26l77, (A26b) 

SEIiy) = 796 - (32 607/ + 1157/^ - I29OO77 - 20707/7/) /(I + OMy) , (A26c) 

dEl{y) = 1500 - 6907/ + 46007/ , (A26d) 



(A25a) 
(A25b) 
(A25c) 
(A25d) 



both in MeV units. The normahzation factors are ff — 1 and 



0.28 - 0.0921/ + 0.23y2 - 0.1b/ + 0.017/ , 
0.196 - 0.487y + 2.02r] + O.USyr] , 
0.0298 + 0.0577 - 0.213/ + 0.23/ . 



(A27a) 
(A27b) 
(A27c) 



There is no big difference between i/^ and Ue, because these Gaussians account for the same 
resonance modes. 
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